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Abstract 

We consider the statistical mechanics of a general relativistic one-dimensional self-gravitating 
system. The system consists of iV-particles coupled to lineal gravity and can be considered as 
a model of N relativistically interacting sheets of uniform mass. The partition function and 
one-particle distitrubion functions are computed to leading order in 1 /c where c is the speed 
of light; as c —* oo results for the non-relativistic one-dimensional self-gravitating system 
are recovered. We find that relativistic effects generally cause both position and momen- 
tum distribution functions to become more sharply peaked, and that the temperature of a 
relativistic gas is smaller than its non-relativistic counterpart at the same fixed energy. We 
consider the large- iV limit of our results and compare this to the non-relativistic case. 
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1 Introduction 



One-dimensional systems of N particles mutually interacting through gravitational forces have 
been of interest in astrophysics for more than three decades. While used primarily as prototypes 
for the behaviour of gravity in higher dimensions, one-dimensional self-gravitating systems (OGS's) 
also approximate the behaviour of some physical systems in 3 spatial dimensions. These include 
the dynamics of stars in a direction orthogonal to the plane of a highly flattened galaxy and 
the collisions of flat parallel domain walls moving in directions perpendicular to their surfaces. 
Furthermore, very long-lived core-halo structures in the OGS phase space are known to exist, 
reminiscent of structures observed in globular clusters, in which a dense massive core in near 
equilibrium is surrounded by a halo of stars with high kinetic energy that interact only weakly 
with the core 

The statistical properties of the OGS are particularly intriguing. Despite extensive study, 
many unanswered questions remain. For example it is not clear if the OGS can attain a true 
equilibrium state from arbitrary initial conditions. Its ergodic and equipartition properties are 
still not well understood. This is primarily because the particle interactions of the OGS (as 
with any self-gravitating system) are attractive and cumulatively long range, in strong contrast 
to typical thermodynamic systems for which such interactions are repulsive and short range. For 
the OGS the macroscopic dynamics does not decouple from the microscopic dynamics, and the 
usual thermodynamic analysis does not apply. 

However there are some established features of the OGS. Rybicki derived in closed form the 
single-particle distribution function in both the canonical and microcanonical ensembles. In the 
large N limit these distribution functions reduce to the isothermal solution of the Vlasov equation. 

All studies to date have neglected relativistic effects. This limitation is understandable since 
no relativistic N particle Hamiltonian was available for analysis. However this situation changed 
recently when a prescription for obtaining the Hamiltonian for a relativistic one- dimensional self- 
gravitating system (ROGS) was given by Mann and Ohta 0. This Hamiltonian can be rigorously 
derived from a generally covariant system coupling relativistic gravity in one spatial dimension 
(i.e. a 1+1 dimensional theory of gravity [Q) to N point particles. In the non-relativistic limit, the 
Hamiltonian reduces to that of the OGS. Although not available in closed form, the Hamiltonian 
can be obtained as a series expansion in inverse powers of the speed of light c to arbitary order. 

We consider in this paper the one-particle distribution function for the ROGS. Our work here is 
a natural extension of previous work on the iV-body problem in relativistic gravity. In three spatial 
dimensions an exact solution to this problem is known for pure Newtonian gravity (and a series 
solutions has been constructed for arbitrary N). In the general theory of relativity dissipation 
of energy in the form of gravitational radiation has obstructed progress toward obtaining exact 
solutions to the iV-body problem even when N = 2. However for the ROGS an exact solution 
to the 2-body problem was recently obtained ||, and generalizations including a cosmological 
constant and/or charge subsequently followed || [7|, |J. These solutions include both an explicit 
expression for the proper separtion of the two bodies as a function of time and an explicit expression 
for the Hamiltonian for the 2-body ROGS as a function of the proper separation and the centre- 
of-inertia momentum of the bodies. 

Encouraged by these results, we here begin a first attempt to understand the basic features of 
the iV-body ROGS. We shall recapitulate the canonical formalism used in ref. || to derive the 
Hamiltonian for the iV-body ROGS. We then compute the partition function and canonical dis- 
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tribution functions. Using an integral transform we then calculate the microcanonical distribution 
functions. All results are in closed form to leading order in 1/c. We consider the limit of large 
N and compare the ROGS and the OGS. We close with a few remarks. Lengthy intermediate 
calculations are confined to appendices. 

2 Canonical iV-particle Hamiltonian of the ROGS 



The OGS Hamiltonian is 



H = H T" 2- + 7TG Y.Y. m a m b \z a - Z b 



a b 
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= T,V^ + 2 ^ G T, m ^b\z a -Z b \ (1) 

a Zm a- a>b 

where the summation is over all N particles, located at positions z a along the spatial axis. The 
potential term straightforwardly follows upon solving the Newtonian equation 

VV = 47rGp (2) 

in one spatial dimension, where p = J2 a rn a 8{x — z a ) is the mass density of the N point particles. 
Our task in this section is to find a prescription for obtaining a relativistic generalization of ([!]) . 

The Hamiltonian for the ROGS that we use is that of a (1 + 1) dimensional theory (a lineal 
gravity theory) that models (3 + 1) dimensional general relativity in that it sets the Ricci scalar 
equal to the trace of the stress-energy of prescribed matter fields and sources. Hence matter 
governs the evolution of spacetime curvature which reciprocally governs the evolution of matter 
[|J . We refer to this theory as R = T theory. Apart from being able to model a number of textbook 
scenarios in general relativity |J, it has the attractive feature of having a consistent Newtonian 
limit [[|. This limit, essential for our purposes, is problematic in a generic (1 + l)-dimensional 
theory of gravity theory JT(| . 

Since the Einstein action is a topological invariant in (1 + 1) dimensions, a scalar (dilaton) field 
must be included in the action pl| . Its coupling to the curvature is chosen so that only the trace 
of the stress energy of matter (N point particles here) is set equal to the Ricci scalar. This action 
will form the basis for the ROGS we consider. Upon canonical reduction of the action ||, the 
ROGS Hamiltonian is given in terms of a spatial integral of the second derivative of the dilaton 
field, which is a function of the coordinates and momenta of the particles and is determined from 
the constraint equations. 

The action integral for the gravitational field coupled to iV point particles is 

1 , .... r_ 1. 



(3) 
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where \I/ is the dilaton field, g^ v and g are the metric and its determinant, R is the Ricci scalar, 
and r a is the proper time of a-th particle whose mass is m a , with k = 8nG/c 4 . We use V M to 
denote the covariant derivative associated with g^ v . 

The field equations derived from the action (|3|) are 

#-<f V M V^ = 0, (4) 
- V M ^V^ - i^V A ^V A ^ + V A V A ^ - V M V,^ = kT^ + -g^A , (5) 
(I f ; l , ,^4 A 



where 



Tfu, = J2 m a[ dTa—^g^gvp Za f a S 2 (x - z a (r a )) , (7) 

a y 9 (M~ a CLT a 



is the stress-energy due to the point masses. Eq.(j5]) guarantees the conservation of T^ u . Inserting 
the trace of Eq. (|5]) into Eq. (f|) yields 

R-A = kT^. (8) 

Eqs. (H), © and (||) form a closed sytem of equations for gravity and matter. 

In order to obtain the Hamiltonian in canonical form, we first decompose the scalar curvature 
in terms of the extrinsic curvature K via 

^T~gR = -2d (^K) + 2d 1 [(N 1 K - d^/^/j] , (9) 

where the metric is 

ds 2 = -N 2 dt 2 + 7 (^dx + ^-dtj , ( 10 ) 

with K = {2N Ql )-\2d 1 N 1 - 7 - 1 iV 1 «9 1 7 - d ol ), so that 7 = 9n,N = {-g m )- 1/2 and N x = g 10 . 
Rewriting the action (||) in first-order form yields 

1 = Jd 2 x ^ Pa zJ(x - z a (t)) + 7T 7 + 11* + N R° + A^i? 1 j (11) 

where it and II are the respective conjugate momenta to 7 and *. Here 



a V 7 

i? 1 = -7T- -IM' + 2n' + Y — S(x-z a (t)) , (12) 

7 7 a 7 

with the symbols ( ' ) and ( ' ) denoting <9 and d±, respectively. 
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Variation of the action (|ll|) yields the set of equations 

* + N {^^-^U + -±— (*?-£ ^—S(x-z a (t)) 

2 ^7 8*^77 a 2 ^M + ml 

+ * l-- 2 W + - + T,- 2 ^- *.(*))) + N o^-*' + N[- = , (13) 

7-iV (2K v ^y 7 7r-2/€ v ^n) + iVi^--2iV{ = 0, (14) 

7 

i?° = , (15) 

R l = , (16) 

ri + ^(-I^n + — !— w + -^aq = o , (17) 

7 2k^ ft^ 

* + iV (2KV77r) - ^i(-*0 = 0, (18) 

7 



dN p* 2 iVo ^7 Wp a 
Pa + -5— a/ — +m 1 -To -~ 

oz a V 7 2 M + m 2 7 ^ a ^ 7 



4 - iVp 7 + = . (20) 
/£| + m 2 7 



All metric components (N , Ni, 7) in equations (|19|) and (|20"1) are evaluated at the point 
where 

df _ df(x) 



dz a dx 

The quantities Nq and N\ are Lagrange multipliers which yield the constraint equations (|15|) and 
fll6|). The above set of equations can be proved to be equivalent to the set of equations (|), (Q) 
and © 1. 

An examination of the generator of space and time transformations [§, 0] indicates that we 
find that we can consistently choose the coordinate conditions 

7 = 1 and II = (21) 

upon which the action fly]) reduces to 



1 = Jd 2 x ^£p a z a 8{x - z a ) - n\ , (22) 



after elimination of the constraints, where 



H = I dxH = — / A* (23) 



is the Hamiltonian for the ROGS. 

The field \I/ is no longer arbitrary, but is instead a function of z a and p a that is determined by 
solving the constraints which are now 

1 N , 

- -(^') 2 + «V +^^ + m 2 a S{x -z a ) = 0, (24) 

^ a=l 

AT 

2tt' + $> a 5(z - z a ) = . (25) 

a=l 

once the coordinate conditions fl2T|) are imposed. Equation (|24D is an energy-balance equation 
which states that the energy of the particles plus the (negative) gravitational energy must vanish. 
Equation (|25|) states that the total momentum of the gravitational field and the particles must 
vanish. The consistency of this canonical reduction was proved in ref. |5| by showing that the 
canonical equations of motion derived from the reduced Hamiltonian (|23| ) are identical with the 
equations ([19|) and (|2~0D . 

Hence the Hamiltonian ( p3| ) is a function only of the coordinates and momenta of the N 
particles in the system. We turn next to its evaluation. 



3 Computation of the ROGS Hamiltonian 



Although the constraint equations are straightforward to solve in the regions between the particles, 
the matching conditions of these solutions at the juncture of the particles are quite non-trivial. 
For the 2-body ROGS their enforcement yields an equation which determines the Hamiltonian in 
terms of the remaining degrees of freedom of the system. While this procedure holds in principle 
for the iV-body ROGS, we have not found a tractable means of obtaining an analogous determining 
equation for the Hamiltonian. 

However it is possible to straightforwardly and rigourously construct approximation schemes 
for computing the ROGS Hamiltonian for iV particles. For example, the post-linear approximation 
is an expansion of the Hamiltonian in powers of the gravitational coupling k, obtained by writing 



= ^(D + ^(2) + 

X = 



x (0) + kx {1) + 



(26) 
(27) 



wherex is defined by x' = tt- Insertion of these expansions into eqs. ( p4| , p5"D yields 



H(2) = Y.\/P 2 a + m2 a+ %Y.Y.[\/pl + m2 a"JPb + m b-PaPb)\rab\ 
a ° a b 



+ ? Z)S ( ^Pl + m lPb ~ Pay/Pb + m b ) 7 ~ab 



w a b 
1 //^ 2 



+4I4J <£vrf+™2 



J2Pb \ r ab\ + € J2\/Pb+ m b r ab 
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E Pb M + C J2\Pb + m b r ab 



E M + 



m 2 c | rod 



+ e E^ 



+ EE \/ Pi + m l\J Pb + m I Kb\ -£Pa\jpl + mlr a b 



-EE 



a 



X 



•c ■ m c 



Pi + mlp b \r ab \ - e p a p b r ab 



\ r bc\ + e E^6c 



E^ |r fc | + e E VPc + ™c r *>. 



(28) 







upon insertion into ([23|), where r ab = z a — z b is the relative separation between particles a and b. It 
can be shown that the solutions to eqs. (|2^ , |25|) must satisfy the boundary condition \]> 2 — 4k 2 x 2 
in the regions | x \»\ z a | in order for the Hamiltonian to be finite 

The ^-expansion is appropriate for describing relativistic fast-motion of the particles and can 
be carried out to any desired order. However to compare the ROGS from R = T theory with 
the OGS, we turn to the post-Newtonian expansion, which is an expansion of the Hamiltonian in 
powers of c . Since both p 2 /m 2 and a/k are of the order of c~ 2 all terms up to the order of c~ 4 



are included in the post-linear Hamiltonian ([28|). The post-Newtonian Hamiltonian to this order 
is therefore [Hi 



H 



N N p 2 N 

m aC 2 + E 7T 2 - + 2nG E m a m b \r a b\ 
a=l a=l Zm a a>b 



N 



Pa , KG 



N N 



2 E o 3 2 E E 771,1 \ r ab\ 



2nG " 

Z^PaPb \r a b\ 

a=l^"°a ^ a=lb=l "° ^ a>b 

fnG\ 2 N N N 
+ ( ) EEE m " m ^ m c [|r ab | |r ac | -r ah r ac ] H 

^ C ' a=l b=l c=l 



(29) 



where the explicit powers of c have been restored. 

The first term in fl2D|) is the total rest energy of the particles, and the second two terms 
are the OGS Hamiltonian ([I]). The remaining terms are all relativistic corrections to the OGS 
to order c -2 . The first of these corrections is a special-relativistic one, whereas the remaining 
corrections are due to relativistic gravity in one spatial dimension. Note that gravity not only 
modifies the potential to a quadratic form, but also includes couplings between particle momenta 
and their positions. These features - modifications of the distance behaviour of the potential and 
position-momentum couplings - are fully analogous to those in general relativity in three spatial 
dimensions. 

We next find the equations of motion for the position of the ath particle. This follows straight- 
forwardly from Hamilton's principle. We have 



dH 

dp a 
Pa_ 

m n 



Pi 
2ml 



k N n k N 

TL m i \ r ab\ + -Tl^Pb \ r ab\ 



(30) 
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and 



Pa 



dH 

dz a 



N 1 r N / p 2 2 \ N 

-2-KGy^m a m b e ab - — ttGV \m a — + m b — e ab - 2-nGy"p a p h e c 



JV JV 



+2(7rG) ^5Z m " m fe m c t^fe kac| +£ab\rbc 



6=1 c=l 



(31) 



where e ab 



1 ^ a > 2 6 
-1 z a < z b 



We can solve (|30|) for p a : 



p a = m a i a + — 



N 



m a -^ + 2ttG^ ((m b m a z a - m a m b z b ) \r ab \) 

1 6=1 



+ 



(32) 



and then insert this into the equation (^T|) for p 



1 

m a z a + — 



3z 2 N 

m a z a -^- - 2nGj2 ((m b m a z a - m a m b z b ) \r ab \ + (m b m a z a - m a m b z b ) e ab (z a - z b )) 

1 6=1 



JV 



-2nG ^a^b 8 



ab 



6=1 



JV 

irG ^ ( m a 



JV 



m b 



6=1 



■2 | -2 

Z b i Z a 



e ab - 2-kG Y m a m b z a z b e ab 

6=1 



JV JV 



+2 (ttG) m a m b m c [e ab \r ac \ + e ab \r bc \ - r ab ] 

6=1 c=l 



(33) 



which simplifies to 



m n z. 



a -a 



N JV 

2txG ^ m a m b e ab + — ^ (m a m b [?>z 2 a + [z b - i a ] 2 } e c 



6=1 



6=1 



^2nG\ 



2 JV JV 



C 
TTG 

c 



) J2J2 m am b m c (e ac - e bc )\r ab \ 

6=1 c=l 

JV JV 

Y m a m b m c [e ab \r ac \ + e ab \r bc \ - r ab ] 



(34) 



6=1 c=l 



upon an iterative substitution of z a in powers of c _1 . These equations of motion reduce to those 
of the OGS in the limit c — > oo , and may be shown to be equivalent to the geodesic equations to 
this order ||. 

We wish to investigate the instrinsic structure of the system described by the Hamiltonian 
fl2Tjp . However because of the translation invariance of the system, two phase-space degrees of 
freedom are redundant, and so must be factored out; otherwise certain average properties such as 
density would be uniform throughout space. 

Using equation(|3T|), it is straightforward to show that 



JV 

£p- = o 

a=l 



(35) 
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and (using also fl30|) ) that the Hamiltonian is time- independent. This means that we can perform 
the phase-space integration subject to the constraint 



p = 



N 

where p = p, 



'a 



(36) 



0=1 



since we can choose a frame of reference in which the centre of inertia is constant. 

Removing the redundant position degree-of-freedom is somewhat more delicate. Although the 
system is invariant under the translation z a — > z a + z, eq. (|32|) implies 



which cannot be written as a total time derivative. Physically, the centre of inertia is the 
relativistically well-defined concept, whereas the centre of mass is not. However we can deal with 
this problem by inserting a factor of unity in all phase-space averages in the form 



where z = YSa=i m a z a , with M = J2a=i m a- If the L-dependence of any integral trivially factors 
out (or can be removed by a shift of variable in the integrand), then we regard the remaining 
quantity as the physically relevant one to describe the system. 

For a canonical ensemble, all phase-space averages are carried out with a weighting function 
exp [—/3H], where fcgT = f3~ l is the temperature multiplied by Boltzmann's constant ks- For 
the microcanonical ensemble, an additional constraint of fixed total energy 



must be included, consistent with the time-independence of the Hamiltonian (|29|) . Since the 
system is in momentum isolation, it is difficult to see how it can be in energy contact with a 
heat bath, and so the physical relevance of the canonical ensemble is somewhat unclear. However 
an evaluation of quantities within the canonical ensemble is instructive in its own right and is a 
necessary preliminary to computing quantities in the more realistic microcanonical ensemble, and 
so we include it in the present discussion. 

Henceforth we set m a = m, so that M = Nm. 

4 The Canonical Ensemble 

We consider in this section the relativistic corrections to the canonical one-particle distribution 
function f^(p, z), which is defined to be the phase-space average of the quantity 




(37) 




H(z a ,p a ) = E 




a 




weighted by exp {—(3H) with the constraint p = 0. Hence 
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Z = — J J dpdz5 (p) J°° dL5{z- L) exp {-(3H) (40) 



where 



is the partition function and where the 2nd line in (p% follows from the indistinguishability of the 
particles. Note that a shift of integration variable 

z' a = z a + L 

renders the partition function in the form 

Z=—J°°dL J J dpdz'5 (p) 5 (z') exp (-/3H) (41) 

where the L-dependence is seen to trivially factor out. It can therefore be dropped (along with the 
prime notation) from further consideration in the evaluation of Z. Similarly, the single-particle 
distribution function becomes 

/f (P, z ) = YN l /_~ dL J J d P dz6 (?) 6 CO ex P N'^Siz-L- z' a ) 5 (p - p a ) (42) 

which is of the form J^dL f^ R (p, z — L). We therefore regard fc R {p,z — L) as the physically 
relevant quantity, where 

f'c R (p, z ) = YnJ J d P dz6 (P) S CO ex P (-^) ^ £ 5 (* - O 5 (P - ^) (43) 

and where the primes will hencefore be dropped. 

Unlike the non-relativistic case, neither the partition function nor f R (p, z) are separable. We 
proceed by first evaluating the partition function. 

We first write the Hamiltonian (|29| ) in the following form 



so that 



H = Mc 2 + H + \h r (44) 

cr 

N 2 N 

H o = £|^ + 27rGm 2 ]T|r af) | (45) 

W ^4 JV JV TV 

Hr = -^2-^ + ^Gj2^2Pb\ r ab\-2nGj2PaPb\r ab \ 

0=1 ° m a=l 6=1 a>6 

/V /V N 

+ (^) 2 E £ £ ™ 3 [M Kd - r ab r a c} (46) 

a=l 6=1 c=l 

exp (-/3if) = c -^ e -^o ^ _ ^ F ^ + W (47) 



which is valid to the order in which we are working. Writing 

N 



& (p) = — [ dk exp 
2n J 



ik E Pa 



n 

0=1 



(45 
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we have 

p -fiMc 2 

z =— 



N 



J dz5 (z) exp | -2-KGpm 2 ^ |r a6 | j J ^ J dp exp (ik ^ p a ~ P X J~ 



(49) 



4.1 The Partition Function 



Consider first the integral 



(50) 



which has integrals that are at most quartic in the momenta. Straightforward Gaussian integration 
yields 



*/ d pexp(«!>-''gl|){ 



PaPb 
4 

I Pa ) 



n{ (3 ) 



(?) 



1 



a = b 



(?) 



2 3(iV-l) 2 



AT2 



(51) 



Hence we obtain 



1 / 27rm 



(iV-l)/2 



P 



N N N 



i - — (ttG) XI m3 tl r ^l l r «cl - ^r ac ] 



a=l 6=1 c=l 



N \ P J 

ttGP (m\ N- 1 J^J^, . 2ttG/? / m \ /3 3 (iV — l) 2 (m 

' 1 2^ 2^ H — | -^7^ J 2^ + 



Nf3 



a>b 



g m 3 c 2 N 2 



2 / AT N 



[3(TrG) 2 """ 3n . . 2irGm". , 3 (N - 1) 2 \ 1 [2nm\ 
1 2 1^1^ [\ r ab\\r a c\-r ab r ac \ — 2jr afe | + 



a=l 6=1 c=l 



a>fe 



0/ Va=l / 

2X 1 ^7rm\ (iV 
8Npmc 2 ) VN \ P ) 



We next consider the integration over the spatial variables. Introducing 

1 N 

Ul = Zl+i - Zl 1 < I < N U N = — Z rn (53) 

we have 



N i N N N-l N-l 

EW = oEEM = E I (Jv-')(^i-fl) = E'(^-')«i (54) 

a>b Z a=l 6=1 1=1 1=1 
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where without loss of generality the particles are ordered in the sequence z\ < z 2 < • • • < Zn, and 
the overall result is then multiplied by N! This gives 

p roc roc roc roc 

/ dz5 (z) F(z) = N! / dzt / dz 2 / dz 3 --- dz N 5 (z) F(z) 

J J—oc J Zl J Z2 Jzn-i 

/co roc roc roc 

duN I du\ I du 2 --- / duN^\5 (un) F(u) 
-oo JO JO Jo 

= N! / rfuF(u) (55) 



provided the function F(z) is symmetric under interchange of any pair of variables, which is the 
case here. The inverse transformation reads 

\ N—l f _ l n>[ 

z n = u N - — D n,i where D n>l = i N _ [ n<[ (56) 

l=n \ — 

and so we have 

N N N I N N 

EEE l\ r ab\ Vac\ ~ rabTac] = 2 E r ba r ac + E Wo 
a=l 6=1 c=l \b>a>c b>a>c 

N-2 N-l 

= 4E E {N-l)(l-k)ku lUk (57) 

k=l l=k+l 



z 



Consequently the partition function is 



N \ (3 J 

ABm 3 
x 1 — 



\(N-l)/2 -oo -oo f oo / N-l 

J J du\ J du2 • • • J duN-i exp f — 2nGf3m 2 n (N — n)u 



N-2 N-l 



2nGm 



N-l 



7T 



G) E E (N-l)(l-k)ku lUk —J2l(N-l) Ul + 



k=l l=k+l 



The three basic integrals in ( |o*%| ) are 

J du exp ( —A ^2 n (N — n) u r 



N-l 



N-l 



n=l 
N-l 



du ^2 k (N — k) Uk exp —A ^ n (N — n) u n 



i=i 



N-2 N-l 



n=l 
N-l 



duY, E (N-l) (l- k) ku lUk exp -A^n(iV- 



n)u r 



k=l l=k+l 



n=l 



yielding 

Z = e 



1 (2ixm\ 

7T) 



(N-l)/2 



( 



-, , 3(iV-l) 2 
~ 8Nf3mc 2 



{2irGfim 2 f~ L [{N-l)\y 



i=i 



3 (N-l) 
8Npmc 2 



1 



\N-i[(N-iyY 

N-l 
X N [(N- l)!] 2 

s~^N-2^N-l (l-k) 
2^k=l l^l=k+l l(N-k) 

\N+1 [(AT - l)!] 2 



(58) 



(59) 
(60) 

(61) 
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4/9 



TV — 1 



(itGfm 3 ^ (l-k) 2nGm 

c 2 {2nG(3m?) N+1 W - l)!] 2 ti l(N - k) c 2 (2KG/3m*f [(N - l)!] 2 



-(3Mc 2 ( 2ttto \ 



(N-l)/2 



N (27rG(3m 2 f- 1 [(N - l)!] 2 



1 - 



1 f (57V + 3) (N — 1) + 87V E^ 1 ESi 



(3mc 2 



8N 



exp 



^Mc 2 - In (/W) - | ( 5JV + 3 )( JV ~ 1 )+ 8JV EL" 1 1 t^t 



SNBmc 



y/N(V2^G/(^y N 1} [(^-1)!] 2 
which is the partition function to lowest relativistic order. The average energy is 

(E) = 

= M C 2 + 1 (jv - 1) /r 1 • 



(62) 



(57V + 3) (TV - 1) + SjVgg Egj^ ■ 
8Mc 2 



(63) 



to the relevant order in c~ 2 . The relativistic correction grows quadratically with N (for fixed 
M = Nm) and is negative. Hence the average energy of a relativistic gravitating system is lower 
than its non-relativistic counterpart at the same temperature. 

4.2 The Single-particle Distribution Function 

Consider next the one-particle distribution function, which is 



N 



JV n=l 



(64) 



where 
flip,*) 



,-/3Mc 2 



C ' 2 f - I N \ 

— J dz8 (z) exp I -2ixGf3m 2 ^ \r a b\ I S(z- z n ) 

/oo /*oo /-ex) poo / 

d^iv y <^ M i y du 2 ---J du^-iS (un) exp ^— 2TcGf3m 2 ^ I (N — l)ui 



x 



JV-1 



X 



5 ( Z - U N + — D n,lUl J 6> cn (P, Z) 



(65) 



e -/3Mc 2 



i-oo poo / N-l \ / ]_ AT-1 \ 

y o dui • • • ^ rfwAf.i exp ^-2rrG(3m 2 J2 C i u ij S y z + J2 D n,iUij cn (p, z) 



in which 



e cn (p, z) = / £ / dp exp ^ f; Pa - p j: ||) - 1# *) 5 (p - p») ( 66 ) 
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and 



where eq. 



Ci = l{N -I) 
j) has been used to express 



Y N-l iv-i j 

Z n = U N + — lu l - E U l = U N - T7 D n,lUl 
iV 1=1 l=n iV 1=1 



N-l 



(N-l) 
-I 



n < I 
n > I 



(67) 



N-l 



(68) 



Evaluation of 9 cn (p, z) is somewhat lengthy, and so we relegate its computation to the appendix. 
We obtain 



cn {p, z) 



: exp 



Nf3p 2 



2iTm 



(N-2)/2 



y/N=T ~" r V 2m (JV - l)J \ (3 J 

f3m^Gmf N N N 



2(3mc 2 



N N N \ 

EEE tl r afe| Kc\ ~ r ah r ac ] 

a=l b=l c=l / 

(3 2 p 4 N (JV 2 - 3JV + 3) 3(3p 2 (JV - 2) 3 (JV - 2) 



(2m) 2 (JV-l) 3 2m(JV-l) 5/2 4(JV-1) 



(69) 



+ 



27rmG\ 
c 2 ) 



1 



JV JV 



JV 



" 0=1 6=1 

or alternatively, in terms of the u variables 

Nf3p 2 \ (2nm\ {N - 2)/2 



2m(N-l) 2 2(iV-l) 



f N 

E 

\c=l 



Ocn(P, U) 



1 



VJV^l 



cxp 



2m(N- 1)J \ (3 J 



C k=l l=k+l 



2(3mc 2 
2nmG 



PYN (JV 2 - 3iV + 3) 3/?p 2 (JV - 2) 3 (JV - 2) 



(2m) 2 (JV - l) d 

7V-1 

. 1=1 



2m (JV - 1) 5/2 4(JV-1) 3/2 
JV 2 /?p 2 JV x /n - x 



(70) 



JV-l 



^m(JV-l) 2 2(JV-l, y v „ , 
Now consider eq. ( |6~5|) , which can be rewritten as 

e-? Mc2 rdk 



Us 



z 



r dk f ( iV " 1 \ 

J — J rfuexp y-ikz — A/3 ^ (Q + iaD nl ) uij 9 cn (p, u 



where 



A = 2^Gm a 



N/3X 



(71) 



(72) 



The integration now involves straightforward integrations over the u variables, after which an 
evaluation of the /c-integral using Jordan's lemma must be performed. This involves some rather 
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tedious manipulations which we describe in appendix ??? The final result is 



y^vrm (JV — 1) 



N-l 

xE 
1 



f3mc 2 
1 



8iV 



fc=i ;=fc+i 



l(N-k) 



(3mc 2 
2 / 



2/3mc 2 ^ (2m) 2 (iV - l) 3 
B? - Af (l - 2iV (/^vrGm 2 ) / |z|) 



/?V (1 + (TV - l)- 1 ) 3(jV-2)/V 3 (iV - 2) 5 



2m(JV-l) 2 4(JV-1) 



A' 



(3mc 2 1^2 (iV — 1) 2m 
1 



iV 2 



(iv-ir 



Cf - jAf (l - 2iV (pnGm 2 ) I \z\) 



(73) 



f3mc 2 



cxp 



+ ^ (l - 2JV (/^vrGm 2 ) Z |*|) 
Integration over p yields the canonical density distribution function 

/oo 
dpfcniP, Z) 
-00 

= (2nGm 2 N!3) exp 



Nf3p 2 



2m (iV — 1) 



- 2nGNpm 2 l \z\ 



1_ j (5iV + 3) (JV-1) ^ ^ (l-fc) 



k=1 i =k+ i 

P + /W 1,8 TV A S| A J 



AT-l 



1 



f3mc 2 



A? - K? 



(l - 2itGm 2 N(3l \z\) 1 exp {-2itGm 2 N(3l \z\) (74) 



whereas integration over z yields the canonical momentum distribution function 

/•OD 

$cn(p) = / dzf cn (p,z) 



(N(3) 



\ 2nm(N - 1) 



cxp 



Npp 2 



2m (JV — 1)_ 



x 1 + 



f3mc 2 



' Nj3 2 p 4 (N 2 — 3JV + 3) 
, 8m 2 (JV - l) 3 



/3 p 2 ( 4jV 2 _ 7jV + g) 57V ^ _ 1) + 3^ 



4m (JV — l) z 



8JV(JV- 1) 



(75) 



From the results in the appendix, it can be shown that 



JV-l 



N-l 



N-l 



E <= E (^-1) E yCf = 

Z=l f Z J=l 4 Z Z=l 4 



AT-l 



JV 



which can be used to show that 

/oo 
dz p c 
-00 



1 and 



/oo 
dp$cn{p) 
-00 



(76) 



(77) 
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We also have the relations 



n-i 1 /V- 1 N ~ x 1 (N- l) 2 

V A? = LSLJl V B» = - { } (78) 

fr[ 1 2 2N-3 ^ 1 2 2N-3 ( J 



the first of which was demonstrated by Rybicki 



5 The Micro canonical Ensemble 



The results for the canonical ensemble obtained in the previous section are for a system of relativis- 
tic gravitating particles coupled to a heat bath which keeps the system at a constant temperature 
T = In such a situation the energy of the system ill-defined, and undergoes fluctuations of 

the order kT. An isolated system, on the other hand, would have its total energy conserved, and 
this is the more realistic astrophysical case. This entails usage of the microcanonical ensemble, 
in which phase space integrations are carried out by constraining the total energy to be E. The 
weighting function e~^ H in the phase space integral is therefore replaced with 5 (E — H). 

Fortunately it is straightforward to compute the relevant microcanonical quantities from the 
canonical ones. Using the same reasoning that led to (f42|), the microcanonical single-particle 
distribution function is 

/£(P> *) = T^n / / d P dz6 (P) 6 00 S(E-H) N- 1 J25(z-z a )8(p- Pa ) (79) 

where 

^ = w\l I dpdz6 (p) 6 ® 5 ^ E ~ H ) ( 8 °) 

Note that Q and Z are related by the Laplace transforms 

Z = / dEe' pH n(E) (81) 
Jo 

n = hl d ^ Ezw (82) 

where in contour C in the latter integral extends from — ioo to -Moo to the right of all singularities. 
Using the general result that 



Mr 1 



it is straightforward to obtain 

(27rm) (iv - 1 » /2 (E - Mc 2 f N ~ 5) ' 2 



n 



VN (vrGm 2 )^ 1 [(N - l)!] 2 r (| (N - 1)) 
x exp 



(E - Mc 2 ) I (5N + 3) (N - 1) + SiVE^i 1 ^ 



(JV-fc) 

Mc 2 1 12(JV - 1) 



54) 
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to leading order in 1/c. 

Similarly using ( |33"D in (]79|) , we have 

Using (|B"2"D and (|73|) we have 

1 



J5) 



^27rm (JV- 1) 



(5JV + 3)(jV-l) + ^ ^ (/ - fc) 



iv-i 



(3mc 2 
2 



11 + 



/3mc 2 | S.Y 

i //?y (i + (at - 1) 



fc=l l=k+ 



1 l(N — k) 



2(3mc 2 \ (2m) 2 (N - l) 3 
5f - A? (l - 2A^ (/37rGm 2 ) / 1 
(3p 2 



3(N-2) (3p 2 3(N -2) 2 
2m (AT — l) 2 + 4(iV- 1) 



A' 



/3mc 2 V,2 (JV — 1) 2m 
1 



N 2 



C C (P, z) 



(3mc 2 
2itGm 2 



(N — I) _ , 
D? + K? (l - 2N (/3vrGm 2 ) I 



Cf - yAf (l - 2N (/?vrGm 2 ) / \z\ 



i6) 



exp 



2m (N - 1) 



2kGN f3m 2 l \z\ 



/ n \ 3/2 r (| (jv - i / 

y/27rm(N-l)) W - Mc 2 ) ) r (§ (N - 2)) 



x exp 

N-l 



(E - Mc 2 ) \ (5AT + 3) (AT - 1) + SJVEjfc 1 £f=*+i jjfcSy 



x 



E 



Mc 2 



N 



3p 2 (A" - 2) 
mc 2 [4m (JV - 1) 2 _ 



+ 



i=i 

2N{7cGm 2 ) \z\ ( , , ( A' 



12(iV- 1) 



4* 



JV 



2^2 



2N 2 p 



/ \+lKi 



N 



(£ - Mc 2 ) 2mc 2 
(E - Mc 2 



(N-l) 

p 4 N (N 2 - 3N + 3) 2p 2 

m 



2m(N - 1) \) _ 

n^3N/2-4 



T 



3JV/2-4 



(fiV-3 



a; 



(2m) 2 (AT - 1)' 
3(N -2f 



2N 3 7rGm 2 



mc* 



+ 1 



(N-l) 
^ * (CI 



T 



3JV/2-5 



V _ l^iV 



(N-l) 



N I x 3iV/2-3 



HN-1) 

as the expression for the relativistic microcanonical partition function, valid to O (1/c 2 ), where 

Np 2 2NnGm 2 



T(p,z) = l 



2m(N -1)(E - Mc 2 ) (E - Mc 2 ) 



I z 



Employing the expressions 



f_ x dy (i - y 



r(7 + i) 
r (7 + l) 



39) 
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f *vv 2 { l -y' 



r(7 + i) 
2r ( 7 + 1) 

3r (7 + 1) 
4r ( 7 + f) 



(90) 
(91) 



the density distribution is 

/oo 
-oo 

/ 2NnGm 2 \ 



(E - Mc 2 ) J 



cxp 



(E - Mc 2 ) \ (5iV + 3) (JV - 1) + SiVE^l 1 Ef^+i 



Mc 2 



2( £ GMj)N f 



Mc 2 



+ f (g :. Mc2) W,))f^ 



Mc 2 



1 



12(JV- 1) 

(T(0,^))f /2 - 7/2 
• AT 2 



A? 



3 (JV 2 - 3JV + 3) 3 (JV - 2) 3iV (JV - 2) 



(iV-1) 



+ 



+ 



4(JV-1) 8(JV-1) 



(N-l) 



+ JV 



+ I -JV5f + 



N 2 (pr - \Af 
{N-l) 



-N[D? + Ki 



' 2NnGm 2 \ 
(E-Mc 2 ) ) 



cxp 



(E - Mc 2 ) \ (5N + 3) (JV - 1) + SJVE^! 1 Ef^i ^Sy 



AT-1 



3JV- 5 



.4 



Mc 2 

N 2 (vrGM 2 ) / |z| 



12(JV- 1) 



Mc 2 V"* 
- 5/2 ^3(JV-l) 2 



(T(0,*)) 



3JV/2-7/2 



Af — JV (Ef - ) - JV (d* + 



The normalization of the density f^o dz p mc (z) = 1 implies 



2 ^|^L + ({E-Mc 2 )\ (A? ({N-l)' 



2NB? 



2N Df 



Mc 2 



I 



4 



exp 



(E - Mc 2 ) J (5JV + 3) (JV - 1) + Ef=*w 



3(JV-1)Z 3(JV-1) / 



j(N-k) 



Mc 2 



12(JV- 1) 



(93) 



which, using (|76|), is easily shown to be satisfied to first order in £ = M 2 ; where the latter 
quantity is the dimensionless fraction of excess energy above the total rest mass. 

The momentum distribution is 

/•CO 

$ mc (p) = / dzf£ c (p,z) 
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N 



\2irm(N-l) (E - Mc 2 ) 
x exp 



\i/2 r (|(iv-i) 



rffiv-2 



(E - Mc 2 ) \ W + 3) (N - 1) + SJVEjfc 1 EjSu -jgy 



Mc 2 



12(iV — 1) 



N-l 



2^ 



A' 



+ 



2Aj 
Imc 2 



+ 



'4* 



3p 2 (JV - 2) 
4m (AT — 1) 2 _ 

A 7 



2(Cf 



1 yjiV 



Imc 2 



N 2 p 2 



m(N — 1) _ 



A(E — Mc 2 
mc 2 (3N - 4) ^l 2 " I (A 7 - 1) 



/v, 



(T(p,0))f/ 2 - 2 



2^ 



.El 

m 



AT 2 



Z 2 mc 2 



(T(p,0))f /2 " 3 



(AT — 1) J J 

N ^p 4 N (N 2 — 3N + 3) 



iV 



(E - Mc 2 ) /mc 2 

4(£-Mc 2 ) / A 
(3N - 4) mc 2 \T 



l 4 AT 



/(iV-1) 
N 



(2m) 2 (N - l) c 
3(N - 2) 2 
8(iV-l) + 

_Dl_+Kl 
I 



3N/2-A 



(T(p,0)) 
/ 

(T(p,0))f /2 - 2 



v 2 r(|(iv-i)) 
r(|iv- 

(JS - Mc 2 ) { (5iV + 3) (AT - 1) + SiVEtY Efix ^ 



\2nm(N -I) (E-Mc 2 ) 
x exp 



Mc 2 
N 



12(N — 1) 



Mc 2 1 4m 



W 2 - 11JV + 6 



(N-iy 



(T(p,0)) 



3AT/2-3 



A 7 



3 (A 7 — 2) 



(E - Mc 2 ) 4Mc 2 



^A 7 (A 72 - 3A 7 + 3) 
(2m) 2 (A 7 -!) 3 



2A 7 (E - Mc 2 ) / (5A 72 - 20N + 12) 



(3A 7 - 4) Mc 2 



(AT-1) 



(T(p,0))f/ 2 " 4 

(t-s) 
htMiHN-s) 



(T(p,0)) 



3AT/2-3 



(T(p,0)) 



3N/2-2 



(94) 



It is straightforward to show that /f^ dp-d mc (p) = 1 from (EO - RTI) to first order in £ = jinr ^ ■ 



6 The Large TV Limit 



For statisical systems (such as those of interest in stellar dynamics), the large A" limit is of 
considerable physical interest. This is the limit in which the total energy E and total mass M = 
ATm are fixed. In the non-relativistic case, the single-particle distribution function approaches 
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the isothermal solution of the Vlasov equations in the large N limit ||. However the relativistic 
case is somewhat more subtle, since the expressions we have obtained are valid only if the speed 
of light is sufficiently large relative to other quantities of the same dimension, and as N becomes 
large we must ensure that this approximation remains valid. 

In order to investigate the large N limit it is necessary to rewrite all quantities in terms of E, 
M, and N. As in ref. H, we adopt the dimensionless variables 



V 



V 



where 



L 



mV 

2 (E - Mc 2 ) 2(c 2 



V 



2 _ 4 (E - Mc 2 ) _ 4(c 2 



(95) 



(96) 



3ttGM 2 3itGM 3M 3 

are the characterisitic length and velocity scales of the system. The scaled distributions functions 
are correspondingly defined 



so that 



■§*(r]) = mV-Q(mVr]) 
f'* R (ri,£) = mVLf' R (mVr},L£) 



i 



(97) 



(98) 



Consider first the partition function (|62|), which can be rewritten as 

-[3Mc 2 



cxp 



Z 



(5N+3)(N-l)+8N^- 1 1 J2l- 1 +1 -V^- 



i(JV-fc) 



(3Mc 2 



N 



3(JV-l)/2 



— / , \ (N-l) 

N(V2^G/c 3 ) [(JV-1)!] 



(3Mc 2 



(99) 



The approximation (K7|) is valid provided 



^N-l (j-fc) 



R , (jjg + 3) (iV - 1) + SiVE^j 1 Ef =fc+1 & 
1 8Mc 2 [ } 

which sets an upper bound on the thermal energy kT = (3~ x of the system for a given value of N. 
For fixed Mc 2 the relativistic corrections are valid only as the temperature becomes vanishingly 
small in the limit of large N. The exponential approximation is slightly better than the polynomial 
one because of the positivity of the partition function. 



Consider next the average energy in the canonical case as given by eq. (|5"3"D, which we rewrite 

as 

(C) = g _ _J_ _ f (»T + 3) (JV - 1) + SN E ^ ffl. a im' (101) 
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Figure 1: Maximum value of the average relativistic energy as a function of N. 

where C s l£^!l, as before . When the thenna. e„e rgy *T = of the 8y ste,„ is su ffl eie„t. y 
small relative to its rest energy Mc 2 , the expression for the average energy (0 does not differ 
much from the non-relativistic value given by its first term. As the thermal energy grows (i.e. as 
/3 decreases) the value of (0 increases more slowly than its non-relativistic counterpart, reaching 
a maximum when 

R (jjg + 3) (N - 1) + 8iVEg ESi lg§) - a 

P ~ 6 (N - 1) Mc 2 = ^ (1 ° 2) 

after which the average energy decreases with decreasing j3, becoming negative when (3 = |/3m- 
As N becomes large, (3m — > f| — ^f~) At /3 = /?m the average energy has half the value of 

its non-relativistic counterpart. In figure(l] we plot the maximum value of (0 as a function of N. 
The curve asymptotes to the constant value of {0 = .573940872 as iV — > oo. 

Of course the relativistic expansion ([0]) breaks down well before f3 reaches this point. In 
figure we plot the average energy (0 as a function of kT for N = 10. The relativistic case 
is clearly distinguishable from its non-relativistic counterpart once kT '/Mc 2 > .02. However the 
upper bound on the thermal energy is kT /Mc 2 < .0117 for N = 10. In ^| we plot the average 
energy over the allowed range of kT illustrating that the distinction between the two cases is about 
8% at most. For iV = 1000 the maximum difference between the two cases is less than one part 
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Figure 2: Average Energy (C)as a function of kT for N = 10 for the non-relativistic and relativistic 
cases. Axes are in units of Mc 2 . 
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;ure 3: Average energy (() for N = 10 over the allowed range of kT. 
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in a thousand over the allowed range of kT. 



In the canonical case we take the energy E to be the fixed total average energy as given by eq. 
3|). Solving this equation for the inverse temperature (3 yields 



18 (JV-l)- (5JV + 3) + ^ EtY EKi 



l(N-l) 



c 



12 (E - Mc 2 ) 



Qjv 
Mc 2 C 



— C 

°jv / 



(103) 



where 



JV-l JV-l 



- 2 (N-1) 



(3 - k) 



(104) 



are defined for convenience. The limit (|100|) implies that 

c / 



\a N a N j 



Cn 



4 9 / s 



(105) 



where the latter limit holds for large JV. Since the exponential forms a better approximation than 



the polynomial one, the value of £ max is probably a bit larger than what is given in eq. ( |105|) 
although it is not clear how much. We obtain 



M 



N 



VLf R {mV^LO = ^\/f (C<?) 3/2 f ,R (mV V ,LZ) 
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(106) 
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Figure 4: The non-relativistic canonical density function for various values of N. 

which is valid only to first order in (. 
The canonical density (|T4T ) becomes 




exp 



-2 (1 - l/N) 1 



(107) 



We plot in fig. f| the non-relativistic canonical single-particle density function p*(£,C — 0) for 
various values of N, recovering the results of Rybicki0. With the exception of iV = 2, the central 
density grows with increasing iV and the distribution becomes slightly more sharply peaked. It 
can be shown that as iV — > oo the canonical density p*(£,C — 0) — > \ sech 2 ^, and the single 
particle distribution function approaches the isothermal solution of the Vlasov equation H. 



In figs. ^| - H we plot the relativistic canonical function p*(£; C) f° r differing values of N. As 
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Figure 5: The canonical density function for N = 3 for various values of the relativistic paramter 
(. The non-relativistic red curve has ( = 0, and the blue, green and black curves are respectively 
( = 0.1, ( = 0.3, and ( = 0.5. Here eq. (|105|) yields C max ~ 0.43 as the limit for which the 
relativistic expression is valid. 

is clear from these figures, relativistic effects significantly enhance the central density by as much 
as 30% depending on the magnitude of (. Even for ( = 0.3, the central density is larger than its 
value of 1/2 in the non-relativistic large N limit. The falloff of the relativistic density functions 
is also more rapid than in the non-relativistic case. 

Unfortunately there is no closed-form expression for the terms Kf and Df and so it is 
not possible to evaluate an explicit expression for either p* (£) or p* mc (£) in the large N limit. 
Instead these quantities must be computed using symbolic algebra; for iV > 20 this involve the 
factorization of thousands of terms, and computer memory limitations make this a prohibitive 
task. However the large N behaviour should not be too different from the N = 20 case, at least 
for small values of £■ Figures | and |10| plot the density for different iV at two different values of 

C- 

The canonical momentum distribution function flTSP is straightforwardly evaluated to be 



/4c c z 

$*cn(v) = mV$ cn (aVr]) = md ^-$ cn (aVr]) 
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Figure 6: The canonical density function for N = 5 where the red blue, green and black curves 
respectively correspond to ( = 0, 0.1, 0.2, 0.3, and ( mSuX ~ 0.29. 
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Figure 7: The canonical density function for N = 10 where the red blue, green and black curves 
respectively correspond to ( = 0, 0.1, 0.2, 0.3, and ( mSuX ~ 0.15. 
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Figure 8: The canonical density function for N = 20 where the red blue, green and black curves 
respectively correspond to £ = 0, 0.1, 0.2, 0.3, and C max ~ 0.15. 
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Figure 9: The canonical density for different values of N at ( = 0.1 
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Figure 10: The canonical density for different values of N at ( = 0.3 
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Figure 11: The canonical momentum density as a function of rj for N = 3. On the left hand side is 
the behaviour of 9* (77; () and on the right hand side is its behaviour relative to the non-relativistic 
density 9* (77; 0). The curves are respectively red, blue, green and black for ( = 0, 0.1, 0.2, 0.3. 
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where each of ( |1071 , |108| ) are also valid to first order in (. Here we can easily take the large N limit, 
which is 

PM - + ^(^^ - V - 9 - 2tt 2 (2r? 2 - l))e^ 2 

to leading order in (. 



(109) 



We plot in figs. [11] - [T^ the behaviour of the canonical momentum density as a function of 
the rescaled momentum r\ for differing values of N. The central momentum density increases with 
increasing (, and falls off more rapidly than in the non-relativistic case. 

However for 77 > 2, the momentum density grows relative to its non-relativistic counterpart, 
overtaking this value for large enough 77 . The relative growth is exponential, although the overall 
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Figure 12: The canonical momentum density as a function of rj for N = 30. Notation is as in fig 
TSL 
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Figure 13: The canonical momentum density as a function of r] for iV = oo. Notation is as in fig 
TIL 



momentum density is exponentially damped for any £. As N increases, the differences between the 
non-relativistic and relativistic cases become less pronounced, although the basic features remain 
the same even in the limit that N — ► oo. 
The microcanonical results are 
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where each of ( p.lO|Jlll|J112| ) are valid to first order in (. As N — > oo we have 



if 



^ i(l + ^(18^ 4 - 81^ + 27 - 2tt 2 (2^ 2 - l))e"^ 



(113) 



Note that this differs from the canonical momentum density unless ( = 0. 

The micro canonical density function Pm C (£,i C) f° r differing values of N is plotted in figs p^-4T7|. 
For N = 3, the microcanonical density is uniform until £ = 9/8, after which it falls linearly to 
zero. For £ < 9/8, at most two particles can contribute to the density; in this region relativistic 
effects enhance their contribution. However for £ > 9/8, at most one particle can contribute, and 
relativistic effects suppress its contribution until £ = 9/4, after which the density vanishes. For 
larger N, as in the canonical case, relativistic effects significantly enhance the central density by 
as much as 30%, depending on the size of (. Their falloff is more rapid, and for sufficiently large 
£ the non-relativistic density distribution is larger. As N increases, the density becomes more 
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Figure 14: The microcanonical density function for N = 3 for various values of the relativistic 
paramter (. The non-relativistic red curve has ( = 0, and the blue, green and black curves are 
respectively ( = 0.1, ( = 0.2, and ( = 0.3. 
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Figure 15: The microcanonical density function for N = 7 for various values of the relativistic 
paramter (. Notation is as in fig. [14]. 
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Figure 16: The microcanonical density function for N = 10 for various values of the relativistic 
paramter (. Notation is as in fig. [14]. 
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Figure 17: The microcanonical density function for N = 20 for various values of the relativistic 
paramter (. Notation is as in fig. [14]. 
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Figure 18: The microcanonical density for different values of iV at ( = 0.1 



sharply peaked and the contrasts between the non-relativistic and relativistic cases become less 
pronounced. 



In figs. [18], [19| we plot the microcanonical density distribution for increasing values of N and 
fixed £• 

We plot the microcanonical momentum distributions in figs. |20H22[ The results are qualita- 
tively similar to the canonical case, although the actual functional forms differ. 

For small N, the relativistic approximation breaks down even for ( as small as 0.3, and the 
momentum distribution function goes negative, as shown in fig.^. However for larger iV the 
momentum distribution is positive for all ( < 0.3, for example as in fig. The relativistic 

densities are more sharply peaked and for sufficiently large momentum parameter r] are larger 
than their non-relativistic counterparts. 



7 Closing Remarks 

We have carried out the first analysis of the statistical behaviour of a ROGS to leading order in 
1/c. The qualitative behaviour of the ROGS as compared to its non-relativistic OGS counterpart 
PI is clear. At a given energy, the ROGS temperature is smaller than the OGS temperature; 
relativistic effects cool the gas down. The one-particle distribution functions become more sharply 
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Figure 19: The microcanonical density for different values of N at £ 



40 



1 



Figure 20: The microcanonical momentum density as a function of r] for N = 3. On the left 
hand side is the behaviour of Q* mc (77; () and on the right hand side is its behaviour relative to 
the non-relativistic density 9* (r/;0). The curves are respectively red, blue, green and black for 
C = 0,0.1,0.2,0.3. 
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peaked in each case with increasing N. For a given N, the ROGS density functions become more 
sharply peaked as the relativistic parameter ( increases. For both canonical and microcanonical 
distribution functions, POGS > ^ROGS ^ or sufficiently large position parameter £. However for 
the momentum densities, this is not true. Although $OGS > ^ROGS ^ or intermediate values 
of the momentum parameter rj, once rj becomes large enough this inequality is reversed. This 
behaviour is presumably due to quadratic character of the ROGS potential relative to its linear 
OGS counterpart in the former case, and from the p 4 corrections in the Hamiltonian fl46|) in the 
latter situation. 

This work can be extended in several directions. It would be straightforward to extend these 
results to the charged and cosmological systems considered in refs. [7J [J to see what effects these 
impose on the distribution functions. Extensions to unequal masses would also be interesting, 
although considerably more difficult. It would be very interesting to go beyond leading order in 
1/c to investigate non-perturbative effects of the ROGS. 

Further understanding the ROGS will undoubtedly require numerical experiments for various 
values of N. The equations of motion yield quartic (as opposed to quadratic) time-dependence 
of the position variables (to leading order in 1/c), and so can be straightforwardly (although 
somewhat tediously) integrated to investigate its equilibrium and equipartion properties. Work 
on this is in progress. 
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8 Appendix 

8.1 The Gaussian integrations 

The following Gaussian integrations were used to obtain (|5T|) 




(114) 
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of which the relevant part is 
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where (1571) was used. 



8.3 The Ai part 

Now we have 
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8.3.1 The A 2 part 

Now we must do 
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8.3.2 The A 3 part 

Now we must do 
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8.3.3 The full expression for 9 cn (p,z) 
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8.4 An evaluation of f cn (p,z) 

Now consider eq. (??), which can be rewritten as 
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This function has single poles at a — in, where n is an integer taking its values between —N and 
N, except for N = 0. Its residues are 
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The remaining integrals are somewhat more difficult. We obtain 
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its values between —N and N, except for N = 0. Writing 

N-l -\ 

^(N -ia) - (1 - ia) = 



ia 



we have the last line of ( |146| ) proportional to 
r(l -ia) r(l +ia) 



r (N - ia) T(N + ia) 

r (1 - ia) T (1 + ia) 
r (TV — ia) r (TV + ia) 



7V-1 



iV - 1+ia 

s=l 
AT-1 

iV - 3 - 2a 2 



s=£n 



ia 
1 



s + ia 



s 2 + a z 



2n 2 



a 2 + n, 2 



where the term in curly brackets contains the double-pole at a = in. Hence 



1 N f dk f ( 7V_1 \ N ~ 1 

-T7 E / 7T~ / duexp -ikz - X(3J2 ( c i + iaD m) uA ^ s(N - s 

iV n=l J Z7r J \ 1=1 ) s=l 



Us 



(147) 



NpX r da 



2tt 



(A/?) 



N 



-iaNf3\z 



r(l- ia) T(l + ia) 
r (TV - ia) T (AT + ia) 



iV 



JV-l 

3 - 2a 2 £ 



s 2 + a z 



2n 2 



a 2 + n 2 



The residue from the double pole at a = in is 



. d 

dx 



2n 2 x(-l) n exp(N[3\xz) 
r (N - x) T (N + x) (x + n) 



—in (—l) n exp(NpXnz) 



[1 + 2N(3\nz + 2n (^ (N - n) - * (N + n))] 



2T (N - n) V (N + n 
provided z < 0, and so the total residue at a = in from (|147 ) is 



(148) 



in (— l) 71 exp(A^/5Anz) 

r (iv - n) r (N + n) 

—i exp(N(3\nz) 



K N-l 

N- b - + 2n^ 



s 2 — n 2 



[T(N)Y 



B 



N 



A% (1 + Nf3\nz 



+ n (* (N + n) - m (N - n)) - (1 + iV/3An,2) 

(149) 
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giving finally 



1 N r dk f ( N ~ 1 \ N ~ x 

E / 7T / ^uexp -ikz (°i + iaD m) ui) J2 S (N 

iV n=l J Zn J \ 1=1 / s=l 



S) U s 



2A i,iG '" 2 X [B? -A?(l- 2N ( hCnr) I \z\)] exp ( 2.Y ( hCnr) I fl)o) 



(27iG(3m 2 f [T{N)] 



where 



B 



N 



n(-l) n+1 [T(N)Y 
T(N -n)T(N + n) 



N-l 



n--+2u 2 y: 



s 2 — n 2 



+ n(V (N + n) - * (N - n)) 



(151) 



Using the relation 



N 

E 

c=l 



n-l 7V-1 

X>,+ £(jv-z) Uj 

Z=l Z=n 



(152) 



we find 
1 

N 
1 

iV 



AT-l 



/ rfuexp -A/3 ^ (Q + iaD n i) u, ^ \r 9 

n=l J \ 1=1 

N (n-l N-l \ 



n=l \s=l 
1 



2 \N 



AT {2-nGpm 2 ) 
1 



E 

71=1 



n-l 



E 



j (N — s — ia) 



N-l 

E 



i 



s + ia) 



/iV-l 



n 



Try— — -7v E{[^( 1 - Ar + ^)- ,I '(^- Ar + ^) + ^( Ar + ^)-^(^ + ^)] 
A {2nG(3m z ) n=1 



1 



x 



r (N - n + 1 - ia) T (n + ia) 



(n-l)\(N-n)\ T(N-ia) T (N + ia) 



1 



N 



AT , na 2 nv E {[* - *«) + * ( N + ia ) ~ * - « + 1 - *«) " * (« + 
N (27iGpm 2 ) n= i 



1 



x 



r (N - n + 1 - ia) T (n + ia) 



1 



N 

V 

N (2irGf3m 2 ) N £=\ 



(n-l)\(N-n)\ T(N-ia) T (N + ia) J J 

1 r (N - n + 1 - ia) F (n + ia) 



(N -ia) + ^(N + ia)] 



d 



1 



(n-l)\(N-n)\ T(N-ia) T (N + ia) t 
T (AT — n + 1 — ia + <r) T (n + ia + cr) 



(27rG/5m 2 ) 



2 \^ 



(AT — ia) + ^(N + ia)] 



da \(n - 1)\(N - n)\ T (N — ia) 

r(l-ia) r(l + ia) 



r (AT + ia) 



J (T=0 



r (JV - ia) r (AT + ia) 
1 d (T(N +l + 2a)T(l-ia + a)T(l + ia + ay 



T(N+l)da\ 
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T (2 + 2a) r (N — ia) T(N + ia) 



J <7=Q 



r(l-za) r(l + za) 



(2nG{3m 2 T (N - ia) r (iV + ia) 

x[V(N- ia) + * (JV + ia) - * (1 - ia) - * (1 + ia) - 2* (JV + 1) + 2^ (2)] 



r(l-ia) r(l + ia) 



(27rG/5m 2 ) T (N — ia) T (N + ia) [f^ s- ia 



N-l 

E 



i Ar - 1 i 
+ E 



JV-1 



r(l-ia) r(l + za) 



(27rG(3m 2 ) N T (N — ia) T (N + 



ia 



s=l 
JV-1 



s + ia 



2E 



1 



AT-1 -i 

j s — ia j s + ia 



s 

N-l 



-2E 



1 s + !. 
1 



s=l 



s + 1 



2n 



a 2 + ri 2 



(153) 



where the curly bracket contains the double-pole at a = in, and all other poles are in the same 
locations as before. 
So we obtain 

1 N r dk r ( N ~ 1 \ ( n ~ 1 N ~ l \ 

jy E J ^ J dnexp y~ ikz ~ X(3 J] (° l + iaDnl ^ Ul ) \£< lU[ + J-' <kN ~ ^7 



N(3\ 



(2nGpm 

~N-1 

E 



da 



2 \ N I 2ix 



e 

N-l 



■iaN/SXz 



r(l-ia) T(l + ia) 
T(N- ia) T(N + ia) 



(154) 



N-l 



+ E 



, s — ia "T s + ia 

s=f=n s^tn 



The total residue at a = in from 



2E 



is 



_T\ n p Nji\nz 



in (—1) e 



\s + l 



N-l 

2E 



+ 



2n 



a 2 + n 2 



s 2 — n 2 



T (N — n)T (N + n) 

+ (^(N + n)-^(N-n)) 
—i exp(Nf3\nz 



N-l 



2E 



1 



8=1 

1 + N(3Xnz 



s + 1 2n 



n 



[T(N)f 



C 



N 



n 



A" (1 + NpXnz) 



(155) 



provided z < (otherwise it vanishes) and so 



1 N f dk f ( N ~ l \ / n_1 N ' 1 

-t-j E / 7T- / duexp -ikz - X/3J2 ( C i + iaD m) uA (J2 lu i + J2( N ~ l ) u i 

iV n =l J Z7T J \ 1=1 / \Z=1 «=n 



2Nj3irGm 2 



N-l 



(27rG(3m 2 f [T(N)] 2 l=1 



Cf - -Af (l - 2iV (>rGm 2 ) / exp (-2JV (>rGm 2 ) Z 



(156) 



where 



n_ n(-l) n+1 {T(N)r 
T (N — n)T (N + n) 



C 



N-l 

2E 



s 2 — n 2 



N-l 



2 — 7-7/ + 777/ + (* {N + n) - * (JV - n)) 



^ s + 1 2n 



(157) 



Finally, we consider the expression 

N „ / N-l 



1 

iV 



AT-1 7V-1 



I rfuexp j -\(3 ( C i + i( x D ni) «i ) E E (N-t)(t-s 

71=1 V Z=l / S=l t=S+l 



SU s U t 
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N 



i N N-l N-l 

EE E (N-t)(t-s)sll S!t 

n=l s=l t=s+l 

1 l N ~ l (t-s)s (N-t) 



(27rG/?m 2 ) N+1 N ^ ^ ^ C s + ia,D ns C t + iaD nt f- = \ Q + iaD nl 

Interchanging the order of the sums gives 

"-i (t-s)s (N-t) 1 
^! C s + ia,D ns C t + ia£> nt C, + iaD nl 

N-l N-l N / u „\ _ / a 7" j.\ /t-i/aA\2 



(158) 



(f - s) s (AT — i) 



(r(AT)) 2 r(A^-n + 1 -ia) r (n + ia) 



SJ^xSv^ + ia^n-Ct + ia^nt^-^KTV-n)! V {N - ia) V{N + ia) i 



T (N + ia) T (N — ia) ti t=T + i 



E 



s) s 



r (N - n + 1 - ia) T (n + ia) * 



+ E 

n=s+l 

+ E 

n=t+l 



^{(N - s)(t + ia)(s + ia) V (N - n + 1) 
r (JV - n + 1 - ia) T (n + ia) \ 



T(n) 



v (iV-s-ia)(t + ia) r(JV-n + l) r (n) / 

(N-t)(t-s) V (N -n + l -ia)T (n + ia)' 

t{N -s-ia){N -t-ia) T [N - n + 1) F(n) 



(159) 



When summed over n, we obtain 



N N-l N-l 

EE E 



(t-s)s (N- t) 



N-l 

n 



v->JV-l 



2k 



N ^ t i^ 1 C s + ic«D„ s C t + iaL> nt ^ C, + iaL> n * fig' 1 (« 2 + ^) 2 



(160) 



which has double poles and single-pole residues at a = ±in for every nonzero value of n < N. 
The coefficients of the polynomial in the numerator are calculable, but we have not found any 
closed-form expression for them. The table below contains results for values up to A^ = 10. 



AT \r^N—l 2k 
iV l^k=l a k a 

3 2(a 2 -l)(a 2 -2) 

4 4(180 - 109a 2 + 10a 4 + 11a 6 ) 

5 8(36576 - 3820a 2 - 75a 4 + 1590a 6 + 149a 8 ) 

6 48(5263200 + 1132124a 2 + 162455a 4 + 170877a 6 + 25445a 8 + 899a 10 ) 

7 288(1455926400 + 635262768a 2 + 123441248a 4 + 43494899a 6 + 6982689a 8 + 399833a 10 + 7163a 12 ) 
1152(1067349830400 + 638760596688a 2 + 149678407480a 4 + 34350170141a 6 + 5098185940a 8 

+351491854a 10 + 10505180a 12 + 110317a 14 ) 
41472(143590977331200 + 103719257351424a 2 + 27508922447056a 4 + 5307728339928a 6 
+709200726957a 8 + 52255261172a 10 + 1966842318a 12 + 35272476a 14 + 237469a 16 ) 
414720(108807366682828800 + 89224919703007488a 2 + 25843463092699920a 4 
10 ; +4662765631167688a 6 + 573351004521465a 8 + 42828568506933a 10 + 1816855681610a 12 

+42140421618a 14 + 493478205a 16 + 2266273a 18 ) 
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Hence we obtain 



N 



N I N-l \ N-l N-l 

E / duexp -\(3 E (Ci + iaDrd) «i £ £ "*)(*" 

n=l \ Z=l / s=l t=s+l 



2 x 2AT/57rG'm 2 



AT-l 



(27rG/3m2 ) N+1 [T(N)} 2 l=1 



E [A N + K i { l - 2N (pnGm 2 ) I \z\)] exp (-2N (finGm 2 ) I 



where the coefficients Df and are determined by the residues given above. These are given 
in the next two tables. 



D? 
= 1 
= 2 
= 3 
= 4 
= 5 
= 6 
= 7 
= 8 
= 9 



N = 3 N = 4 N 



n 

12 



4 
209 



125 

13* 



87 



1500 



45 
169 

$78 



72' 



123480 



N = 6 

31 

144 
5599 

3087 
12781 



mi 



333396 
1005251 

3333960 



N = 7 

2117 

_ 4 Mi 

1344 
1468 



2141 



m 



113400 
22515781 

30187080 
162182479 

1207483200 



N = 8 

2917 

3600 
30347 



2430 



189; 



8085285 
31621187 

25874640 
8830831883 

23686076700 
108510757181 

1989630442800 



N = 9 

24869 

22050 
941 

1125 
836537 



31 



If 



1197900 
440479867 

263178630 
2430681577 

3582153575 
7762820713 

46056260250 
53405900137 

2579150574000 



N = 10 

113931 

78400 
601906 



o11b 



«8 



40 



146210350 
4721124703 

2292578288 
10937730746 

10746460725 
22140852323 

65502236800 
9952938913257 

»» 
20274186832099200 



K? 
= 1 
= 2 
= 3 
= 4 
= 5 
= 6 
= 7 
= 8 
= 9 



N = 3 N = A N = 5 

_1 _1 _3 

,5 



24 



2 

_ J_ 

i 5 
50 



10 
27 

245 
257 

3920 



N = 6 

_2 

w 

784 
829 

10584 
2761 

10584 



N = 7 

5 



224 
11 

1512 
361 

3 ffl 
3659040 



N = 8 

3 



10 

27 
37 

540 
853 

_m_ 

2760615 
248029 

77297220 



AT = 9 
_7 

19§ 

450 
381 

3025 
551 

21780 
44893 

9 » 



1431430 
1797 



234 



32207175 
1075199 

1030629600 



N = 10 

_4 

281? 

605 
1297 

7260 
1153 

204490 
24587 

572572 
66629 

2147145 
274399 

22902880 
1812235 

661893232 
6514549 

19856796960 



Summarizing: 



i « f dk r 



2N(3nGm 1 



N-l 



{2itG!3m 



- \ N ~' [T(N)] 2 l , 
dk 



E A? exp (-2N ^TiGm 2 ) I 



N 



N r dk r f N ~ l 



-jfe2-A/3^^ =1 1 (C(+iaD ni )Mi) 



(162) 
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2N(3nGm 2 



N-l 



(2^G^m 2 f [T(N)} 2 £ 



Y, [B? - A? (l - 2N ((3nGm 2 ) I \z\)] exp (-2N ((3nGm 2 ) I \z\) (163) 



N 



N r Ah r /n-l N-l \ t , 

E / f- du £ < s + £ (TV - S )/0 e(— (Q + ^H) 

n=l" / Z7r J \s=l s=n / 



2A^7rGW 



JV-1 



(2-kG^ 2 ) iV [r(iV)] 2 ^ L 



E ^ " 7^ (l - 2iV (/?7rGm 2 ) Z 



exp (-2N ((3irGm z ) I \z\) (164) 



N 



dk 



'N-l N-l 



,E ^ du £ £ (iV - t) (t - s ) ^ ) el 

iv n=l J Z/l J \ s=l t=s+l 



2N(3nGm 2 



N-l 



]T [D? + Kj* (l - 2N (f3nGm 2 ) Z \z\)] exp (-2JV (/^Gm 2 ) / |s|)(165) 



(2vGPm 2 ) N+1 [T(N)} 2 ti 
The final expression for the one-particle distribution function is 

fcn(p,z) 



VN^G/c^liN-m 2 ^ 



{N - 2)/2 2Nf3nGm 2 (f3mc 2 ) ^ 



yW^T (2irG(3' 



POT 



x exp 

JV-1 

xE 

1=1 

+ 



1 j3(iV-l) 2 (TV - 1) (A 2 (N - 1) + A 3 ) 



(3mc 2 



8N 



-Ai 



^ ^ (Z - fc) 



f w / Ai / W(l + (*-!) 
\ 1 { 2Prnc 2 [ (2m) 2 (N - l) 3 
1 //(A 3 -A 2 )/3p 2 (A 3 + (iV-2)A 2 )' 



3 (N — 2) /3p 2 3(N - 2f 



f3mc 2 \ m {N - I) 2 
2 (2A 3 + (AT - 2)A 2 ) / 1 



(3mc 2 



(N-l) 
Npp 2 
2 (N-l) ~~ 2m 



(JV-l)' 



2m(A^-l) 2 4(iV-l) 
5f — Af (l — 2N (f3iiGm 2 ) I \z 



Cf - -A* (l - 2N (pnGm 2 ) I \z\) 



Df + (l - 2N (/^vrGm 2 ) Z |*|) 



exp 



Npp 2 



(3mc 2 
(2vGm 2 ) (N(3f 2 
y^vrm (N-l) 

1 (3(iV-l) 2 (TV - 1) (A 2 (N - 1) + A 3 ) 



2m (N- 1) 



- 2nGN(3m 2 l \z\ 



x exp 



f3mc 2 



8N 



TV 



JV-1 

*E 

Z=l 



Ax //3V(l + (^-l) d J 3(iV-2)/V 3(JV-2) 



^ ^ (l-k) 
4 fri iti l(N -k) 



2/3mc 2 ^ (2m) 2 (AT — l) 3 2m(iV-l) 2 4 (AT — 1) 

/3mc 2 V m(N- I) 2 (N-l) J L ' 'V 7 ' 

Cf - 1^ (l - 2iV (^vrGm 2 ) Z |z|) 



2 (2A 3 + (N — 2)A 2 ) 



f3mc 2 



2(N-l) 2m L(iV-l) 2 . 
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A i 



(3mc 2 



D? + K? (l - 2N (/JvrGm 2 ) I 



exp 



Nf3p 2 



2m (N — 1) 



2%GNpm 2 l \z\ 



(166) 



The canonical density distribution function is given by integration of f cn {p, z) over p 



Pc 0) 



dpfcn(p, z) 
{2nGm 2 ){N(3) 3/2 \2itm{N-l) 



^2nm (N — 1) V #0 



x exp 



1 J 3 (iV — 1) ; 



(3mc 2 



Ai- 



(iV-l)(A 2 (iV-l) + A 3 ) 



iV-l 

xE 



IN N 
Ax /3(l + (iV-l) 3 ) 



4 h Jki l ( N - k ) 



2j3mc 2 



+ 



3(JV-2) 3(N-2f 



AN 2 (N-l) 2N(N-1) A(N-l) 



/?mc 2 \N(N - 1) 
2 (2A 3 + (jV - 2)A 2 ) vrGm 

(3mc 2 



2(N - 1) 



2(JV- 1) 



Cf - yAf (l - 2iV (/3vrGm 2 ) / |; 



(3mc 2 
[2-nGm 2 N(3 



D? + (l - 2N (/JvrGm 2 ) Z |z|)] 1 exp {-2nGNpm 2 l \z\) 



x exp 

N-l 



f3mc 2 



3 ( -V - l) 2 ^ (N-l) (A 2 (N - 1) + A 3 ) A4 ^ ^ (/-*) 



x 



1 



f3mc 2 



8N 



N 



N 

(A 2 (jV-l) + A 3 ) ' 
N 



k=l l=k+l 



l(N-k) 



B? - \ 4 D? 



f3mc 2 



• (a 2 (N-i) + \ 3 y 

N 



A? - \ A K? 



(l - 2N (/3vrGm 2 ) / 1 exp (-2nGm 2 Nf3l |; 



(167) 



It is straightforward to show that the coefficients Af and Bf obey the sum rules 

N-i 1 AT _ 1 

E 4* - 

i=l 
iV-l 



E^ 



2 2iV-3 

1 (N — l) 2 

2 2iV-3 



(168) 
(169) 



where eq. ( |168j ) was previously derived in the c — > oo limit |2j . Since /f^ cfz p c (z) — 1, we must 
have 



iV-l 



2E K+ 



1 /3Ai(iV-l) 2 



/5mc 2 



iV 



.4 



N 2 / (A 2 (jV-l) + A 3 ) ' 
2 I TV 
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= exp 



f3mc 2 



3{N-iy 
8N 



Ai 



(iV-l)(A 2 (jV-l) + A 3 ) 
N 



JV-l N-l 

a 4 E E 

k=i i=k+i 



or alternatively, to the relevant order in c, 

N-l i 1 

V -A N - - 

i=i 1 z 

N-l i i 

E 5f = Un-i) 

i=i L z 



JV-1 



1 



N-l N-l 



- k =i i=k + i HN-k) 
each of which can be straightforwardly verified. We also can show 



iV-1 



E -cf = cf = (iV " 1} 
i=i i 



N 



E^ 

i=i 



N 



7T 

12 



l(N-k) 



(170) 



(171) 
(172) 
(173) 



(174) 
(175) 
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